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Introduction 


Elastomer  and  rubber  components  are  found  on  many  U.S.  Army  vehicles,  including  tires, 
bushings,  and  track  pads.  An  improved  understanding  of  the  failure  process  in  elastomers 
would  assist  materials  scientists  in  the  design  of  new  and  stronger  rubberlike  materials  and 
would  aid  structural  engineers  in  the  design  of  rubber  components  with  improved  service  life. 
The  ability  to  model  simple  fracture  behavior  imder  large  strains  in  a  viscoelastic,  hyperelastic, 
and  incompressible  material  is  an  important  step  towards  modeling  crack  growth  in  elastomers. 

When  viscoelastic  materials  undergo  a  constant  strain,  the  applied  stress  diminishes  with 
time.  Although  viscoelasticity  is  known  to  affect  crack  tip  behavior  in  hyperelastic  and  incom¬ 
pressible  materials,  its  contribution  has  not  been  quantised.  In  this  investigation,  a  model  for 
nonlinear  viscoelastic  material  behavior  in  a  Mode  I  crack  was  examined.  A  finite  element  anal¬ 
ysis  of  a  Mode  I  Single  Edge  Notched  (SEN)  elastomer  specimen  under  relaxation  conditions 
was  performed  and  compared  to  experimental  data.  Changes  in  the  stress  field  with  time  were 
examined.  The  tearing  energy  was  also  studied  as  a  function  of  time. 


Tearing  Energy 


The  tearing  energy,  T,  was  first  defined  by  Rivlin  and  Thomas  [1]  as  the  energy  released 
per  unit  area  as  a  crack  of  length  L  advances  by  6L  Mathematically,  the  tearing  energy  is 
represented  as 


where  U  is  the  total  strain  energy  stored  through  elastic  deformation  and  b  is  the  specimen 
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tliickness.  When  the  relecise  of  stored  strain  energy  is  greater  than  the  energy  required  to 
generate  new  surface  area,  the  crack  will  propagate.  The  tearing  energy  is  equivalent  to  the 
J-integral[2]  for  nonlinear  elastic  materials.  For  a  crack  advancing  in  the  xi  direction,  refer  to 
Figure  1, 

J  =  jf  [Wm  -  Uui,^\dS  (2) 

where  W  represents  the  strain  energy  density,  ni  is  the  component  of  the  unit  normal  in  the 
xi -direction,  ii  is  the  nominal  traction  vector,  and  Ui  is  the  displacement  vector.  The  J-integral 
is  path  independent  when  evaluated  along  any  suitable  contour  F  in  the  reference  configuration 
which  encompasses  the  crack  tip,  while  dS  is  an  element  of  arc  length  along  F.  For  hyperelastic 
materials,  the  J-integral  is  globally  path  independent  and  equals  the  energy  release  rate.  For 
viscoelastic  materials,  catre  must  be  taken  in  the  interpretation  of  the  J-integral.  The  global 
path  independence  of  J  applies  only  when  the  stress  is  a  single  valued  function  of  strain.  For  a 
Mode  I  crack  under  relaxation  conditions,  this  condition  is  met.  When  the  applied  stresses  are 
infinitesimal,  the  J-integral  is  a  function  of  the  stress  intensity  factor,  K/,  or 


3fi 


(3) 


under  plane  stress  conditions,  where  fi  is  the  infinitesimal  sheau'  modulus. 


Experimental  Procedures 

The  material  used  in  this  study  was  a  carbon  black  filled,  highly  saturated  nitrile  rubber 
(HNBR)  cured  with  peroxide  and  a  zinc  dimethacrylate  co-curative.  It  is  based  on  a  formulation 
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Table  1:  MATERIAL  FORMULATION 


Ingredient  phr 


HNBR  (Zeptol  2020“)  100.0 

N-121  Carbon  Black  20.0 

Zinc  Oxide  2.0 

Agerite  Resin  0.5 

Zinc  Dimethacrylate  coagent®  30.0 

Dicumyl  Peroxide  (Dicup  R)*^  1.3 


“  Zeon  Chemicals,  ACN  36%,  unsaturation  10% 


R.  T.  Vanderbilt,  Polymerized  1,2  -  dihydro  -  2,2,4-  trimethylquinoline 
®  Sartomer,  SR  -  634 
**  Hercules 


developed  for  use  in  tank  track  pads  [3]  and  was  supplied  by  Zeon  Chemicals.  Its  formulation 
is  given  in  Table  1. 

Experimental  tests  were  performed  using  an  Instron  model  4505  screw  type  machine  with 
an  Instron  model  elastomer  extensometer  attachment.  The  use  of  an  extensometer  to  mea¬ 
sure  strain,  rather  than  crosshead  displacement,  eliminates  measurement  error  from  slippage. 
Pneumatic/hydraulic  grips  were  used  with  pressures  from  7  to  10  MPa.  Load  and  extension 
data  were  electronically  captured.  Strip  specimens,  measuring  15.24  by  2.54  cm,  were  die  cut 
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from  standard  tensile  sheets.  The  specimen  was  subjected  to  conditioning  between  40%  and 
200%  tensile  strain  at  a  strain  rate  of  2%/sec  for  10  cycles.  The  conditioning  was  performed 
to  reduce  the  Mullins  effect.  After  the  specimen  was  unloaded  and  removed  from  the  grips,  it 
was  allowed  to  relax  for  15  minutes.  A  crack  of  length  0.64  cm  was  then  introduced  along  the 
mid-height  of  the  specimen  with  a  razor  blade. 

After  conditioning  and  introducing  the  crack,  the  specimen  was  reinserted  into  the  test 
machine.  The  specimens  were  pulled  at  a  crosshead  rate  of  0.5%/sec  to  50%  strain  and  the 
load  history  was  recorded  for  15  minutes.  In  order  to  accurately  capture  the  peak  stress,  the 
data  was  recorded  at  a  speed  of  10  points/sec  for  the  first  five  minutes.  Thereafter,  the  rate 
at  which  data  was  recorded  was  decreased  to  two  points/sec.  The  specimen  was  then  unloaded 
and  removed  from  the  grips. 


Nonlinear  Viscoelasticity 

The  mathematical  basis  [4,  5]  for  finite  strain  viscoelasticity  will  now  be  reviewed,  focusing 
on  the  Prony  series,  which  is  used  to  model  large  strain  viscoelastic  material  behavior.  The 
implementation  of  the  Prony  series  in  ABAQUS  [6],  the  commercial  finite  element  software 
used  in  this  analysis,  will  also  be  presented.  Consider  a  spatial  or  deformed  coordinate  system, 
y(T),  at  time  T,  and  a  material  coordinate  system,  x,  in  an  isotropic,  homogeneous  and 
incompressible  material.  The  deformation  gradient  is  represented  by  F,  where 

dyij) 

5x 


4 


A  relative  deformation  gradient  at  time  T  relative  to  time  t  cm  be  defined  as 


F(r,t)  = 


MI) 

dyit)  ’ 


(5) 


where 


F(r,t)  =  p(r)F-Hf)  .  (6) 

Under  the  assumptions  of  incompressibility,  the  determinant  of  F  equals  one.  The  left  and  right 
deformation  tensors,  B  and  C,  are  respectively 

B  =  FF^  , 

and  (7) 

C  =  F^F  , 


while  the  strain  invariants,  Jj,  are 


Ii  =  trC  , 


(8) 


For  incompressible  materials,  I3  equals  one.  The  Cauchy  strains,  E,  are  given  in  terms  of  the 


right  deformation  tensor,  C,  as 

E=^|C-I|  ,  (9) 

Here,  I  is  the  identity  matrix. 


The  deviatoric  stress  tensor,  Se,  is  defined  as 


dW.  . 


(10) 
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where  W*  is  the  elastic  strcdn  energy  density.  For  hyperelastic  materials,  can  be  expressed 
as  a  function  of  the  strain  invariants.  Common  examples  include  the  Rivlin  constitutive  model 
and  the  Neo-Hookean  material  law. 

The  Cauchy  stress,  S,  is  represented  by 

S  =  Se  •  (11) 

The  indeterminate  pressure,  p,  is  a  consequence  of  incompressibility  and  is  found  through  the 
applied  boundary  conditions  and  equilibrium.  The  nominal  stress,  tr,  can  be  defined  in  terms 
of  the  Cauchy  stress  as 

<r  =  SF-’’  .  (12) 

The  viscoelastic  Cauchy  stresses  are  now  expressed  as 

/t  ,  _ 

^F-‘(T;i)^r-T(T,()Jr ,  (13) 

In  this  equation,  the  viscoelastic  strain  energy  density,  Wv.  is  a  measure  of  the  energy  dissipated 
by  the  material  and  must  be  measured  from  creep  or  relaxation  experiments.  Consequently,  the 
function  will  decay  or  fade  in  time  and  can  be  described  by  a  monotonically  decreasing 

function.  Although  many  potential  forms  of  are  possible,  a  common  assumption  is 

that  W„  is  the  product  of  separable  functions  of  time  and  strain  throughout  the  relaxation 
spectrum.  A  Prony  series  can  be  used  to  model  the  viscoelastic  response,  as  demonstrated  by 
Quigley  [7],  Johnson  [8],  and  Schapery  [9]  so  that 

dW  ^ 

=  Z)  Prnexp{t  -  T)/t„,  ,  (14) 

m=l 

where  each  time  constant,  t^,  is  paired  with  a  multiplicative  scalar.  Pm- 
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Constitutive  Model 


The  finite  element  analysis  was  performed  with  ABAQUS  [6]  which  assumes  that  Wi,  has  the 
same  form  as  Wg.  The  multiplicative  constants,  Pm,  are  provided  as  nondimensional  fractions, 

^mi  of  a  hyperelastic  energy  function,  W,  so  that 

dW 


where 


and 


Pm  -  9m 


dWe  _  dW 

Qlongterm  ^  j,  > 


SE 


dW„ 

dE 


K 


dW 


=  9kexp[t  -  r /Tk]-^ 


k=l 


It  is  also  assumed  that 


9longterm  "f"  ^  ^  9k  —  1  • 
Jb=l 


(15) 


(16) 


(17) 


(18) 


Therefore,  W  specifies  the  overall  shape  and  magnitude  of  the  combined  elastic  (long  term) 
and  viscoelastic  material  behavior.  The  proportion  of  the  material  behavior  that  is  viscoelastic 
is  governed  by  J^k-i9ki  while  the  proportion  of  long  term  material  behavior  is  determined  by 
9iongterm-  Here,  W  was  arbitrarily  selected  as  the  Rivlin  constitutive  law,  or 

(19) 

where  Ii  and  I2  are  functions  of  the  stretch  ratios,  Aj, 


=  E  E  -  3)’V2  -  3)^’;  C-oo  =  0  , 

i=o  i=o 


Ii  =  Aj  +  A2  +  A|  , 


(20) 


h  =  A?Ai  +  A^Al  +  A?A5  , 
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and  where  Cij  are  time  independent  material  constants.  In  uniaxial  extension,  the  stretch  ratios 
can  be  expressed  eis 


Ai  =  A  , 


(21) 


A2,3  =  1/A^/*  . 


At  an  arbitrary  time,  t,  the  expression  for  nominal  uniaxial  stress  in  a  step  -  strain  relaxation 
test,  based  on  Equations  (10),  (12),  (13),  (16),  and  (17),  becomes 

a(A,t)  =  2(A-l/A^)  •^  +  VA-^  j  '  (22) 

As  t  -♦  0,  substituting  Equation  (18)  into  the  above  expression  for  uniajcicil  stress  yields 

+  )  ,  (23) 

and  stress  strain  behavior  can  be  approximated  from  a  “quick  pull”  test.  Because  o’(A,0) 
is  only  a  function  of  W,  the  material  constants,  Cij,  can  now  be  found.  Previous  studies  [7] 
determined  the  material  constants  for  this  HNBR  elastomer  so  that  the  resulting  constitutive 
model  would  be  stable  in  a  Drucker  sense.  These  time  independent  material  constants  were 
given  the  following  values: 


Coi  =  l.OSMPa  and  Czo  =  0.76  x  IQ-^MPa.  (24) 

A  Prony  series  was  determined  from  relaxation  data[7]  for  this  material.  The  values  of 
the  multiplicative  constants  were  agmn  constrained  to  assure  Drucker  stability  and  are  listed 
in  Table  2.  A  comparison  of  predicted  material  behavior  and  experimental  data  is  shown  in 
Figure  2  for  a  0.5%/sec  single  cycle  hysteresis  loop  to  200%  strain. 
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Table  2:  TIME  DEPENDENT  MATERIAL  CONSTANTS  IN  PRONY  SERIES 


rk 

gk  X  10“* 

0.1778 

5.18 

0.3162 

4.04 

1.0000 

7.19 

1.7783 

2.47 

3.1623 

4.91 

5.6234 

3.35 

10.000 

7.23 

31.623 

7.00 

56.234 

2.88 

177.83 

1.61 

316.23 

8.46 

9 


Mechanical  Crack  Tip  Stress  Field 

There  is  no  theoretical  asymptotic  solution  for  a  Mode  I  crack  in  a  viscoelastic,  hyperelastic 
and  incompressible  material.  Therefore,  for  this  analysis  the  far  field  loads  and  displacements 
from  the  finite  element  analysis  will  be  compared  to  the  experimental  results.  Because  exper¬ 
imental  data  close  to  the  crack  tip  could  not  be  measured,  numerically  determined  crack  tip 
field  quantities  cannot  be  verified. 

Asymptotic  crack  tip  fields  have  been  derived  for  the  mechanical  crack  tip  Mode  I  stress 
field  in  hyperelastic  and  incompressible  materials.  For  completeness,  both  the  linear  elastic 
crack  tip  field  and  the  Mode  I  crack  tip  field  in  a  Neo-Hookean  material  will  be  presented.  At 
the  crack  tip,  a  polar  coordinate  system  was  introduced  such  that 

zi  =  rcosd,  X2  =  rsm0.  (25) 

The  linear  elastic  asymptotic  solution  [10]  is 

^  JiTf  0  /  .  e  .  39\  ,  wjx 

i  5 f ) + ^  ■  (26) 

Kj  0  30  .  0  I  i/jN 

where  Oij  represents  the  nominal  stresses,  and  Kj  is  the  stress  intensity  factor.  This  solution 
would  be  present  in  the  crack  tip  region  only  if  the  far  field  strains  were  infinitesimal  and 
viscoelastic  stresses  were  absent. 

A  Neo-Hookean  material  is  based  on  the  linear,  one  term  form  of  the  Rivlin  law,  Equation(19), 

W  =  C7io(/i  -  3)  .  (27) 
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The  plane  stress  nonlinear  asymptotic  crack  tip  field  for  this  material  was  derived  by  Knowles 
and  Sternberg  [11]  as 


Sn 

S22 

<v 

M  2  -1 

rr  ‘ 

S12 

sin 

2 

yi 

jcosO  , 

0 

2/2 

1/2  •  ^ 
'sm-  , 

(28) 


where  the  Cauchy  stresses,  Sij,  and  the  deformation  field,  yi,  are  referenced  to  the  poleir  coor¬ 
dinate  system  (r,  6)  in  the  undeformed  configuration.  The  load  amplitude  constants,  a  and  &, 
are  functions  of  the  applied  load  and  the  geometry. 


Finite  Element  Model 

A  finite  element  model  of  the  SEN  test  specimen  is  provided  in  Figure  3.  Only  the  top 
half  of  the  specimen  was  modeled  due  to  reflective  symmetry  about  the  Xi  axis.  The  crack 
tip  region  is  modeled  by  a  semicircular  mesh  containing  both  refined  and  coarse  regions.  This 
mesh  design  in  the  crack  tip  region  can  sustmn  large  deformations  at  the  crack  tip  so  that 
nonlinear  stresses  can  be  more  accurately  determined.  At  these  large  deformations,  numerical 
mapping  methods  typically  break  down,  and  the  elements  begin  to  evert,  making  the  nonlinear 
crack  tip  field  difficult  to  resolve.  Consequently,  the  finite  element  mesh  was  composed  of  two 
regions,  a  coarse  mesh  in  the  immediate  vicinity  of  the  crack  tip,  enclosed  by  a  refined  mesh. 
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where  accurate  crack  tip  held  quantities  could  be  found.  The  coarse  mesh  circumvented  the 
numerical  problems  described  above  and  allowed  the  nonlinear  crack  tip  field  to  extend  into  the 
refined  mesh  region.  In  addition,  the  coarse  mesh  region  was  sufficiently  small  so  that  it  did 
not  significantly  influence  crack  tip  field  quantities  in  the  adjacent  refined  mesh  region.  Both 
the  coarse  mesh  and  the  refined  mesh  contained  rings  of  eight-noded  isoparametric  elements. 
Within  a  ring,  all  elements  had  equal  angular  extent  and  the  same  radial  length. 

The  circumferentially  coarse  mesh,  as  shown  in  Figure  4(a),  had  three  rings  of  elements  and 
extended  radially  to  2.54  x  10“^  cm.  The  first  ring  was  constructed  of  three  elements  of  radius 
2.54  X  10“^  cm.  In  each  subsequent  ring,  the  number  of  elements  was  doubled  so  that  the  fourth 
ring  contained  24  elements.  Along  each  circumferential  element  ring,  nodal  displacements  were 
constrained  to  enforce  compatibility. 

The  surrounding  refined  mesh,  partially  shown  in  Figure  4(b),  extended  from  2.54  x  10"^  to 
0.254  cm.  Each  of  the  circumferenticJ  twenty-four  rings  in  the  refined  mesh  were  constructed  of 
24  elements.  Element  radii  were  bieised  such  that,  along  a  radius  extending  from  the  crack  tip, 
they  were  equally  spaced  on  a  logarithmic  scale  from  2.54  X  10"^  to  2.54  x  10“^  cm.  Within 
each  decade  cm  unit  of  crack  length,  to  10“”*,  where  m  ranged  from  -3  to  -1,  were  12 

rings  of  elements.  There  were  21  elements  in  the  coarse  mesh  and  144  elements  in  the  refined 
mesh,  for  a  total  of  165  elements  surrounding  the  crack  tip.  Away  from  the  crack  tip  region, 
the  mesh  was  transitioned  to  a  rectangular  grid.  The  entire  mesh  contained  779  elements  and 
2577  nodes,  with  two  kinematic  degrees  of  freedom  at  each  node  and  one  additional  pressure 
degree  of  freedom  at  each  comer  node. 
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Displacement  boundary  conditions  were  applied  along  the  top  of  the  specimen,  as  shown  in 
Figure  5,  as  a  function  of  time.  The  displacements  uniformly  increased  with  time  for  88  sec.  At 
this  time,  the  stretch  ratios  close  to  the  top  of  the  specimen  equalled  1.5,  while  stresses  should 
be  at  maximum  values  throughout  the  test  specimen.  Prom  88  to  900  sec,  the  displacements 
were  fixed,  allowing  stress  relaxation  to  occur.  The  applied  displacement  history  duplicated  the 
one  applied  during  the  actual  experiment. 

Results 

Comparison  of  finite  element  results  to  experimental  data  verified  that  the  finite  element 
analysis  did  simulate  the  experiment.  Figures  6  and  7  s’how  that  the  experimentally  applied 
stretch  ratios  and  stresses  agreed  closely  with  numerical  data.  The  applied  stretch  ratios  based 
on  numerical  results  were  centroidal  values  obtained  from  rows  of  elements  across  the  top  portion 
of  the  finite  element  mesh.  The  applied  stretch  ratios  increased  steadily  with  the  experimentally 
applied  strmn  rate  to  a  maximum  value  of  1.5  at  88  sec,  then  remained  constant  for  the  rest 
of  the  experiment.  Applied  stresses  were  calculated  from  reaction  forces  across  the  top  half 
the  specimen.  These  stresses  reached  their  peak  value  at  88  sec,  when  the  applied  stretch 
ratio  of  1.5  was  attained,  and  then  steadily  decreased,  mimicking  the  experimentally  observed 
viscoelastic  behavior. 

Energy  changes  in  the  mesh  are  found  in  Figure  8.  The  total  energy  in  the  specimen  is 
the  sum  of  the  recoverable  str<dn  energy  and  the  viscoelastic  energy.  The  recoverable  strain 
energy  was  maximum  at  88  sec  and  then  slowly  decreased,  similar  to  the  applied  stresses.  The 
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viscoelastic  energy  quantity  measured  in  ABAQUS  [6]  represented  the  energy  loss  in  the  finite 
element  model  due  to  viscous  effects  (time  dependent  behavior).  The  viscoelastic  energy  in  the 

test  specimen  increased  steadily  until  the  maximum  applied  stretch  ratio  was  reached.  The 

viscoelastic  energy  continued  to  increase  during  the  relaxation  portion  of  the  experiment,  but 
at  a  slower  rate. 

The  J-integral  was  determined  by  ABAQUS[6]  over  16  paths,  where  each  path  was  a  circle 
encompassing  the  crack  tip.  Global  path  independence  was  maintained  over  the  course  of  the 
entire  analysis,  as  there  was  only  a  maximum  difference  of  3  %  between  the  highest  and  lowest 
values  of  the  J-integral  at  any  given  time.  Average  values  of  the  J-integral  as  a  function  of  time 
are  provided  in  Figure  9.  The  J-integral  increased  until  the  maximum  applied  displacement 
was  attained,  and  then  decrecLsed  as  viscoelastic  energy  losses  occurred,  making  less  energy 
available  for  crack  growth. 

The  deformed  finite  element  mesh  at  time  t  =  21  sec,  as  shown  in  Figure  10,  showed  the 
immediate  development  of  a  parabolic  shape  along  the  crack  surface  and  blunting  at  the  crack 
tip  upon  load  application.  The  blunting  became  more  pronounced  with  load  application  until 
t  =  88  sec,  as  shown  in  Figure  11.  After  t  =  88  sec,  the  deformed  shape  did  not  change  while 
the  specimen  was  relaxing. 

Stresses  in  the  refined  crack  mesh  region  were  studied  as  plots  of  radial  distcince  from  the 
crack  tip,  where  r  refers  to  the  undeformed  element  radius,  and  as  plots  of  the  undeformed 
polar  angle,  6.  Only  centroidal  values  of  stress  were  plotted.  Radial  plots  of  S22,  found  in 
Figure  12,  showed  a  change  in  slope  with  increased  time  and  load  application.  At  t  =  21 
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sec  (A  =  1.1),  the  slope  of  S22  was  constant  and  equalled  1/2.  At  t  =  88  sec  (A  =  1.5),  the 
slope  of  S22  changed  with  distance  from  the  crack  tip.  Close  to  the  crack  tip,  the  slope  was 
approximately  one.  However,  with  increasing  distance  from  the  crack  tip,  the  slope  decreased. 
These  results  suggest  that 


S22  ~  o(r  ^/^)  for  Aoo  ~  1.1 


S22  ~  o(r  ^)  for  Aoo  ~  1-5  . 


(29) 


At  A  ~  1.1,  strains  far  from  the  crack  tip  were  close  to  linear  behavior.  It  is  noted  that  the  linear 
522  stress  component  is  also  of  (refer  to  Equation  (26)).  In  addition,  for  S22  in  the 

nonlinear  plane  stress  asymptotic  crack  tip  field  for  a  Neo-Hookean  material,  see  Equation  (28), 
is  also  of  o(r“^).  These  dominant  stresses  appear  to  be  consistent  with  crack  tip  fields  in  which 
viscoelasticity  effects  are  not  included.  During  relaxation,  when  t  >  88  sec  and  A  =  1.5,  the 
slope  of  S22  was  unchanged  while  the  stresses  gradually  decreased  with  time. 

These  results  confirmed  the  significance  of  viscoelasticity  as  a  fracture  arrest  mechanism 
which  constrsuns  crack  growth.  When  a  crack  grows,  the  newly  created  surface  area  becomes 
traction  free.  As  this  new  surface  unloads,  it  undergoes  viscoelastic  softening,  making  less 
energy  available  for  continued  crack  extension.  These  energy  losses  were  demonstrated  here 
by  the  increase  in  viscoelcistic  energy  with  relaxation,  as  shown  in  Figure  8.  At  the  crack 
tip,  documentation  of  the  sustmned  energy  losses  as  a  function  of  time  was  found  in  both  the 
J-integral,  Figure  9,  and  the  S22  stress  component,  Figure  12. 

Examination  of  circumferential  plots  of  stress  as  a  function  of  the  undeformed  polar  angle. 
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Of  revealed  that  maximuni  values  of  all  three  Cauchy  stress  components  were  found  near  the 
undeformed  crack  flank,  behind  the  crack  tip  (refer  to  Figures  13  to  15).  This  behavior 
was  found  in  a  highly  localized  region  surrounding  the  crack  tip  at  small  radii  for  \  =  1.1. 
The  localized  region  continued  to  expand  in  size  until  the  peak  stretch  ratio,  A  =  1.5,  was 
applied.  The  shape  of  the  circumferential  plots  changed  with  radial  distance  from  the  crack 
tip,  indicating  a  transition  in  the  stress  fleld.  Close  to  the  crack  tip,  the  mechanical  crack  tip 
stress  field  was  found.  With  increasing  distance  from  the  crack  tip,  the  stress  field  changed  to  a 
transitional  stress  state  governed  by  the  presence  of  the  crack  tip,  the  specimen  geometry,  and 
the  applied  load.  This  stress  state  was  shown  at  the  larger  radii  in  Figures  13  to  15. 

Assuming  that  material  failure  occurs  at  sites  of  maximum  stress,  these  plots  indicate  that 
failure  is  more  likely  to  occur  close  to  the  crack  surface,  behind  the  crack  tip  and  suggest  that 
crack  growth  is  likely  to  commence  in  this  region.  Experimental  evidence  [12,  13]  of  crack 
growth  in  plane  stress  SEN  and  other  test  specimens  supports  this  hypothesis.  These  studies 
suggest  that  secondary  crack  growth  occurs  above  and  below  the  crack  tip.  A  plane  strain  finite 
element  analysis[14]  also  predicts  that  secondary  crack  growth  should  initiate  close  to  the  crack 
surface. 


Future  Work 

For  a  better  understanding  of  viscoelastic  material  behavior  in  the  fracture  process  of  elas¬ 
tomers,  additional  experimental  and  cinalytical  resecirch  is  required.  Careful  experimental  stud¬ 
ies  of  crack  propagation  on  cracks  of  varying  lengths  are  needed  to  measure  the  tearing  energy 
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at  the  onset  of  crack  propagation,  zis  well  as  the  direction  and  increment  of  crack  growth.  Fi¬ 
nite  element  analysis  of  these  specimens  would  allow  us  to  monitor  changes  in  stress  fields  and 
energy  surrounding  the  crack  tip  before  and  after  propagation. 
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Applied  Displacements 


Crack  Tip 


Figure  3:  Finite  element  model. 
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Crack  Tip 


Figure  4:  The  finite  element  mesh  in  the  crack  tip  re^on:  (a)  shows  the  coarse  mesh  surrounding 
the  crack  tip;  (b)  shows  the  refined  mesh  region  enclosing  the  coarse  mesh. 
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Figure  5:  Displacement  history. 
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Figure  6:  A  comparison  of  the  applied  stretch  ratios  as  a  function  of  time. 
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Load  (MPa) 


Figure  7:  A  comparison  of  the  applied  load  as  a  function  of  time. 
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Figure  8:  A  history  of  energy  changes  in  the  finite  element  analysis. 
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J'integral  (kPa-m) 


Figure  9:  The  J-iutegral  u  a  function  of  time. 
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Figure  10:  Tlie  deformed  mesh  at  t  =  21  sec. 
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Figure  11:  The  deformed  mesh  at  <  =  88  sec. 
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822/shear  modulus 
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Figure  12:  The  itress  component  S23  u  a  function  of  radial  distance  from  the  crack  tip  at 
9  =  71  degrees. 


30 


s22  /  shear  modulus 


Figure  13:  Circumferential  plots  of  S22  *»  »  fimction  of  polar  angle  at  t  =  88  sec. 
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Figure  14:  Circumferential  ploti  of  Su  as  a  function  of  polar  angle  at  t  =  88  sec. 
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0—0  radius  ■  2.80E-3  cm 
□ — □  radius  *  2.32E-2  cm 
0—0  radius  «  7.32E-2  cm 


Figure  15:  Circumferential  ploti  of  Sn  ai  a  function  of  polar  angle  at  t  =  88  lec 


DISTRIBUTION  LIST 


No.  of 

Copies _  Jo 

1  Office  of  the  Under  Secretary  of  Defense  for  Research  and  Engineering,  The  Pentagon, 
Washington,  DC  20301 

Director,  U.S.  Army  Research  Laboratory,  2800  Powder  Mill  Road,  Adelphi,  MD  20783-1197 
1  ATTN:  AMSRL-OP-SD-TP,  Technical  Publishing  Branch 
1  AMSRL-OP-SD-TA,  Records  Management 

1  AMSRL-OP-SD-TL,  Technical  Library 

Commander,  Defense  Technical  Information  Center,  Cameron  Station,  Building  5, 

5010  Duke  Street,  Alexandria,  VA  23304-6145 

2  ATTN:  DTIC-FDAC 

1  MIA/CINDAS,  Purdue  University,  2595  Yeager  Road,  West  Lafayette,  IN  47905 

Commander,  Army  Research  Office,  P.O.  Box  12211,  Research  Triangle  Park, 

NC  27709-2211 

1  ATTN:  Information  Processing  Office 

Commander,  U.S.  Army  Materiel  Command,  5001  Eisenhower  Avenue,  Alexandria  VA  22333 
1  ATTN:  AMCSCI 


Commander,  U.S.  Army  Materiel  Systems  Analysis  Activity,  Aberdeen  Proving  Ground, 
MD  21005 

1  ATTN:  AMXSY-MP,  H.  Cohen 

Commander,  U.S.  Army  Missile  Command,  Redstone  Arsenal,  AL  35809 
1  ATTN:  AMSMI-RD-CS-R/Doc 

Commander,  U.S.  Army  Armament,  Munitions  and  Chemical  Command,  Dover,  NJ  07801 
1  ATTN:  Technical  Library 

Commander,  U.S.  Army  Natick  Research,  Development  and  Engineering  Center 
Natick,  MA  01760-5010 
1  ATTN:  SATNC-MI,  Technical  Library 

Commander,  U.S.  Army  Satellite  Communications  Agency,  Fort  Monmouth,  NJ  07703 
1  ATTN:  Technical  Document  Center 

Commander,  U.S.  Army  Tank-Automotive  Command,  Warren,  Ml  48397-5000 
1  ATTN:  AMSTA-ZSK 

1  AMSTA-TSL,  Technical  Library 

President,  Airborne,  Electronics  and  Special  Warfare  Board,  Fort  Bragg,  NC  28307 
1  ATTN:  Library 

Director,  U.S.  Army  Research  Laboratory,  Weapons  Technology,  Aberdeen  Proving  Ground 
MD  21005-5066 
1  ATTN:  AMSRL-WT 


34 


No.  of 
Copies 


To 


r 


Commander,  Dugway  Proving  Ground,  UT  84022 
1  ATTN:  Technical  Library,  Technical  Information  Division 

Commander,  U.S.  Army  Research  Laboratory,  2800  Powder  Mill  Road,  Adelphi,  MD  20783 
1  ATTN:  AMSRL-SS 

Director,  Benet  Weapons  Laboratory,  LCWSL,  USA  AMCCOM,  Watervliet,  NY  12189 
1  ATTN:  AMSMC-LCB-TL 

1  AMSMC-LCB-R 

1  AMSMC-LCB-RM 

1  AMSMC-LCB-RP 

Commander,  U.S.  Army  Foreign  Science  and  Technology  Center,  220  7th  Street,  N.E., 
Charlottesville,  VA  22901-5396 

3  ATTN:  AIFRTC,  Applied  Technologies  Branch,  Gerald  Schlesinger 

Commander,  U.S.  Army  Aeromedical  Research  Unit,  P.O.  Box  577,  Fort  Rucker,  AL  36360 
1  ATTN:  Technical  Library 

U.S.  Army  Aviation  Training  Library,  Fort  Rucker,  AL  36360 
1  ATTN:  Building  5906-5907 

Commander,  U.S.  Army  Agency  for  Aviation  Safety,  Fort  Rucker,  AL  3636 
1  ATTN:  Technical  Library 

Commander,  Clarke  Engineer  School  Library,  3202  Nebraska  Ave.,  N.,  Fort  Leonard  Wood, 
MO  65473-5000 
1  ATTN:  Library 

Commander,  U.S.  Army  Engineer  Waterways  Experiment  Station,  P.O.  Box  631,  Vicksburg, 
MS  39180 

1  ATTN:  Research  Center  Library 

Commandant,  U.S.  Army  Quartermaster  School,  Fort  Lee,  VA  23801 
1  ATTN:  Quartermaster  School  Library 

Naval  Research  Laboratory,  Washington,  DC  20375 
1  ATTN:  Code  6384 

Chief  of  Naval  Research,  Arlington,  VA  22217 
1  ATTN:  Code  471 

Commander,  U.S.  Air  Force  Wright  Research  and  Development  Center,  Wright-Patterson 
Air  Force  Base,  OH  45433-6523 
1  ATTN:  WRDC/MLLP,  M.  Forney,  Jr. 

1  WRDC/MLBC,  Mr.  Stanley  Schulman 

U.S.  Department  of  Commerce,  National  Institute  of  Standards  and  Technology,  Gaithersburg 
MD  20899 

1  ATTN:  Stephen  M  Hsu,  Chief,  Ceramics  Division,  Institute  for  Materials  Science 
and  Engineering 


35 


No.  of 
Copies 


To 


1  Committee  on  Marine  Structures,  Marine  Board,  National  Research  Council,  2101  Constitution 
Avenue,  N.W.,  Washington,  DC  2041  8 

1  Materials  Sciences  Corporation,  Suite  250,  500  Office  Center  Drive,  Fort  Washinqton 
PA  19034  ’ 


1  Charles  Stark  Draper  Laboratory,  555  Technology  Square,  Cambridge,  MA  02139 

General  Dynamics,  Convair  Aerospace  Division,  P.O.  Box  748,  Fort  Worth,  TX  76101 
1  ATTN:  Mfg.  Engineering  Technical  Library 

Plastics  Technical  Evaluation  Center,  PLASTEC,  ARDEC,  Bldg.  355N,  Picatinny  Arsenal, 
NJ  07806-5000 
1  ATTN:  Harry  Pebly 

1  Department  of  the  Army,  Aerostructures  Directorate,  MS-266,  U.S.  Army  Aviation  R&T 
Activity  -  AVSCOM,  Langley  Research  Center,  Hampton,  VA  23665-5225 

1  NASA  -  Langley  Research  Center,  Hampton,  VA  23665-5255 

U.S.  Army  Vehicle  Propulsion  Directorate,  NASA  Lewis  Research  Center, 

2100  Brookpark  Road,  Cleveland,  OH  44135-3191 
1  ATTN:  AMSRL-VP 


Director,  Defense  Intelligence  Agency,  Washington,  DC  20340-6053 
1  ATTN:  0DT-5A,  Mr.  Frank  Jaeger 

U.S.  Army  Communications  and  Electronics  Command,  Fort  Monmouth,  NJ  07703 
1  ATTN:  Technical  Library 

U.S.  Army  Research  Laboratory,  Electronic  Power  Sources  Directorate 
Fort  Monmouth,  NJ  07703 

1  ATTN:  Technical  Library 

Director,  U.S.  Army  Research  Laboratory,  Watertown,  MA  02172-0001 

2  ATTN:  AMSRL-OP-WT-IS,  Technical  Library 

1 0  Authors 


36 


